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Abstract. We compute the maximum mass of moder- 
ately rotating strange stars as a function of the strange 
quark mass, of the QCD coupling constant, a c , and of the 
bag constant (vacuum energy density), B, in the MIT bag 
model of quark matter with lowest order quark-gluon in- 
teractions. For a fixed value of B, the maximum stellar 
mass depends only weakly on a c , and is independent of 
this coupling in the limit of massless quarks. However, if 
it is the value of the chemical potential of quark matter 
at zero pressure which is held constant, for example at 
the value corresponding to the stability limit of nucleons 
against conversion to quark matter, the maximum mass of 
the strange star is higher by up to 25% for a c = 0.6, than 
for non-interacting quarks, and this may be relevant in 



the discussion of kHz QPO sources. The maximum mass 



The existence of strange stars would shed light on the 
issue of stability of quark matter. It has been suggested 
that a quark fluid, composed of roughly equal number of 
up, down and s trange quarks, may be the gr ound state 
of bulk matter ( Bodmer, 1971; Witten, 1984 ) and a de- 
tailed descripion of su ch strange matter has been given by 
Far hi and Jaffc (1984 ) . This idea has not been universally 
accepted, and it was argued that for realistic values of 
the QCD coupling constant, the phase transition to quark 
matter would occur at densities too high to be of interest 
(e.g., Bethe, Brown and Cooperstein, 1987). 

Strange stars may be very hard to distinguish from 
neutron stars, particularly i n LMXBs, as they may have 
a crust of normal matter (Alcock et al., 1986]; Haense] 



of a non-rotating strange star could be sufficiently high 
to allow an orbital frequency as low as 1.0 kHz in the 
marginally stable orbit. However, for all a c < 0.6, the 
stellar mass cannot exceed 2.6M Q at any rotational pe- 
riod > 1.6 ms. 
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1. Kilohertz QPOs and the mass of LMXBs 



et al., 1986| ). The crust will contribute little to the mass 
(< 1O _5 M0), but is expected to have a thickness sufficient 
to support nuclear burning, including helium flashes re- 
sponsible for X-ray bursts, thus mimicking neutron stars. 
The heat released in conversion (at the bottom of the 
crust) of nuclei into strange matter is directed into the 
strange matter core ( Miralda-Escude et al., 1990]; Haense] 
and Zdunik, 1990| ) . At any rate, steady release of energy 



The discovery of kHz QPOs ( (Btrohmayer et al., 1996 ; van 
|der Klis et al., 1996 ) in several low-mass X-ray binaries 
(LMXBs) has renewed interest in the maximum mass of 
neutron stars, as its value limits the maximum observable 



orbital frequency (Kluzniak et al., 1990). The maximum 
masses for neutron stars modeled with various equations 



from nuclear conversion is very difficult to distinguish from 
the gravitational binding energy released in steady accre- 
tion. 

However, in young radio pulsars the crust thickness 
does make a difference — the crust in strange stars would 
be far too thin to allow the redistribution of angular 
momentum necessary to explain the glitch phenomenon 
(sudden spin- up of a pulsar), therefore glitching pulsars 
are though to be neutron stars and not quark stars (Al 



par, 1987). It is likely that the coalescence of two strange 



of state (e.o.s.) are well established (Arnett and Bowers, 



1977; Friedman et al., 1986; Cook et al., 1994; Salgado 



stars would lead to the contamination of our Galaxy with 
chunks of strange matter, and would thus preclude the 



formation of young neutron stars (Madsen, 198S; Cald 



et. al., 1994), and have been examined in the context of well and Friedman, 1991). These arguments make unlikely 



kHz QPOs flKluzniak, 1998| ). In principle, LMXBs could 
contain strange stars instead (Cheng and Dai 1996), and it 
has been asked whether the observed values of kHz QPOs 
are compatible with the masses of such stars (Bulik et al., 
1999a,b). 



the presence of strange stars in the population of ordinary 
pulsars (including binary Hulse- Taylor type pulsars) or 
their accreting counterparts, the high-mass X-ray binaries. 
However, the presence of strange stars among th e millisec- 
ond pulsars or LMXBs seems to be allowed ( Kluzniak 
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1994; Cheng and Dai, 1996). In these old systems, the 
stellar mass itself may yield clues as to the nature of the 
compact object. 

We are led to consider the maximum mass of strange 
stars in the expectation that LMXB masses will become 
available in the near future, either through a better un- 
derstanding of accretion phenomena, including the ob- 
served quasi-periodic oscillations (QPOs) in the X-ray 
light curve, or through classical optical determinations of 
binary parameters for newly discovered transient sources. 

The maximum frequency of stable orbital motion is 
attained in the innermost (marginally) stable circular or- 
bit (ISCO) allowed by general relativity, and if its value 
for some X-ray source were that of the observed maxi- 
mum QPO frequency — ranging from 1.0 kHz to 1.2 kHz 
in twelve LMXBs — conclusions could be reached about the 
nature of the compact object. [For a review of QPOs see 
e.g. [van der Klis (1998| ).] 

The question, whether ISCO frequencies as low as 
the maximum observed QPO frequencies can be attained 
outside quark stars, has been answered in the affirma- 
tive for rapidly rotating strange stars, with an equa- 
tion of state (e.o.s.) based on the MIT bag model of 
quark matter — strange stars rotating close to the equato- 
rial mass-shedding limit can have ISCO frequencies below 
1 kHz for masses as low as 1.4M Q ( ^Btergioulas et al., 1999| ). 
However, the same question has not yet been answered for 
moderately rotating strange stars, of periods P ~ 3 ms or 
more, for which such ISCO frequencies would imply larger 
masses: M = 2.2M (1 + 0.75j)(1.0kHz// ma2; ), where 
fmax is the ISCO frequency, and j < 0.3 is the dimension- 



less angular momentum (Kluzniak et al., 1990; Kluzniak. 
19980 

There is a compelling reason to consider stars rotating 
at relatively low rates. The period of the compact star in 



the the transient LMXB SAX J 1808.4-3658 (Wijnand 



and van der Klis, 1998), the first of possibly many such 
transients to be discovered, has been measured to be 2.5 
ms. Whether or not kHz QPOs will be discovered in that 
source, its mass may eventually be determined by optical 
studies of the binary companion. 

It has also been argued that the oscillations seen in 
some X-ray bursters imply a stellar rotational frequency 
of ~ 300 Hz flStrohmayer et al., 19971 ) . Finally, it is not 
yet clear whether strange stars can attain the equatorial 
mass-shedding limit because of the unusually high value 
of T/W > 0.2 calculated for their models — note that for 
Newtonian stars the secular instability to a bar-mode de- 
f ormation sets in alrea dy at 0.1275 < T/W < 0.1375 
( Bonazzola et al., 1996 ). For a discussion of other modes, 
also in general relativity, see, e.g., Gourgoulhon et al. 
(1999), and references therein. 

For all these reasons, we considered exact numerical 
models of strange stars rotating with frequencies up to 
700 Hz, and found that, to within a few per cent, the max- 
imum mass of such strange stars is the same as that of the 
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Fig. 1. The maximum mass of a strange star as a func- 
tion of the strange quark mass, for various values of 
the QCD coupling constant. The stellar mass scales as 
the inverse square root of the bag constant, B, pro- 
vided that the quark mass is scaled with B 1 / 4 . Here, 
B 1 = 58.9MeVnxT 3 . 



static configurations. Therefore, we investigate the max- 
imum mass of static (non- rotating) strange stars, which 
till now has been discussed in the context of QPOs only 
in the simplest case of ideal quark gas in the bag model 



(Bulik et al., 1999a). Here, we use the e.o.s. of interacting, 
massive quarks within the MIT bag model of self-bound 
quark matter. 



2. Strange stars 

The maximum mass of quark stars was first derived 
( Brecher and Caporaso, 1976|) for an unusually low value 
of the bag constant. Wittcn (1984) showed that the maxi- 



mum mass of a static strange star is, within the MIT bag 
model (Farhi and Jaffe 1984), 



M„ 



= IMM^B/Bx)- 1 ^, 



(1) 



where B 1 = 58.9 McV-fm" 3 . 

There seems to have been no systematic investigation 
of the maximum mass of the star as a function of the three 
basic parameters of quark matter in the MIT bag model: 
the mass of the strange quark, m s , the bag constant B, and 
the strength of the QCD coupling constant, a c — g 2 /47r. 
Detailed models of strange stars have been constructed 
and the structure of strange stars, including the mass- 
radius relationship has been discussed extensively in the 



literature (Alcock et al., 1986 


; Hacnscl et al., 1986 


; Glen- 


denning, 1989 


; Fricman and Olinto, 1989; Prakash et al. 


1990), however discussion of stellar parameters tended to 



concentrate on the maximum value of the bag constant, 
setting a lower bound on the maximum mass and an upper 
bound on the rotational frequency. 
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In this section, we discuss only static stellar models, 
constructed in general relativity by solving the TOV equa- 
tions ( Oppcnhcimcr and Volkoff, 1939| ). We neglect the 
crust, whose contribution to the maximum mass is quite 
minor for stars composed mostly of self-bound quark mat- 



ter flHaensel etal., 1986J ). Following Farhi and Jaffe (1984), 
we take m s and a c to be renormalized at q — 313 MeV, 
and as a model for quark matter consider a "bag" of 
positive vacuum energy density, B, filled with quarks (of 
two massless flavors and one massive) having interactions 
through first order in a c . The actual form of the thermo- 



dynamic expressions we use can be found in |Hacnscl et al 
(19860 As both the pressure and the density scale with 
the bag constant, TOV equations imply that the stellar 
mass and radius scale as M oc B~ 1 / 2 1 and R oc £?~ 1//2 



( Witten, 1984 ), provided that the masses of the quarks 
scale as to ex B 1 / 4 . 

In Figure 1, we plot the rescaled value of the maximum 
stellar mass, M(B/ Bi) 1 / 2 , as a function of the rescaled 
strange quark mass, m s (Bi/ B) 1 / 4 , for various values of 
the QCD coupling constant. Note that the highest value 
of M max , the maximum mass of the strange star, is inde- 
pendent of a c , if B is fixed, and is obtained for m s = 0, 
i.e., for massless quarks. But in fact, the value of B is not 
known, and as its lower bound does depend on the value 
of a c , the actual physical bounds on the maximum mass 
of a strange star will depend, through B, on the coupling 
constant. This is discussed in the next section. 



3. The minimum value of B and an upper bound 
to the mass of moderately rotating strange 
stars 

To determine the highest possible value of the maximum 
mass of static strange stars in the MIT bag model, it is 
enough to consider the e.o.s. of an ulrarelativistic Fermi 
gas in a volume with vacuum energy density B > 0. If the 
quarks are massive, the actual maximum mass of the star 
will be somewhat lower — as is evident from Figure 1, at 
fixed values of the other parameters, the maximum mass 
of a strange star decreases with increasing quark mass.f] 

Currently, the actual value of B cannot be reliably de- 
rived from fits to hadronic masses of the quark-model of 
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Fig. 2. The dependence on rotational frequency of the 
maximum mass of a strange star, as derived from se- 
quences of fully relativistic numerical models for the e.o.s. 
of massless, noninteracting quarks (in a bag). In the fig- 
ure, the stellar mass is scaled with the maximum mass 
of static models for this e.o.s., 1.98M . Note, that for all 
periods of stellar rotation in the currently observed range 
(P > 1.6 ms) the maximum mass of rotating strange stars 
differs from the static one by less than 4% (the same result 
holds also for the e.o.s. of strange matter with massive, in- 
teracting quarks [Gondek-Rosinska, et. al 2000]). 

i.e., that the baryonic chemical potential at zero pressure 
of such quark matter satisfies 



> 939.57 MeV. 



(2) 



This condition provides a lower bound on the bag con- 
stant B, and consequently, by eq. (0), an upper bound 
on the mass of static strange stars. We take the up and 
down quarks to be massless. Therefore, eq. @ implies 
B > £>i(l — 2a c /7r). For noninteracting and massless 
quarks (a = 0, m s = 0), this value of the bag constant 
corresponds to a minimum density of strange matter (at- 
tained at zero pressure) of po(0) = 4.20 x 10 14 g/cm 3 = 
ABi/c 2 , with a corresponding maximum mass of non- 
rotating strange stars of M max (0) = 1.98M Q . 

For interacting (but massless) quarks, through lowest 
order in gluon exchange, the equation of state is identical 



nucleoi is. Instead, we must rely on a different argument 
to find B m i n . We require that neutrons do not combine to 
form plasma of deconfined up and down quarks, or equiva- 
lently, that quark matter composed of up and down quarks 
in 1:2 ratio is unstable to emission of neutrons through 
the reaction u + 2d — > n (e.g., Haensel, 1987; Farhi, 1991), 



to that of non-interacting quarks (Chapline and Nauen- 
berg, 1976), p — (p — po)c 2 /3, the only difference being 



1 The dependence of M on m s is easily understood in the 
limit a c — 0. The mass of a static star modeled with ideal 
Fermi gas e.o.s. scales with fermi on mass as M oc m~ 2 and 
attai ns a finite limit at m — * (Oppenheimer and Volkoff. 



193E) 



in that the lower bound on the density at zero pressure, 
following from conditions of neutron stability (eq. |Q ) , is 
decreased with respect to the value for an ideal Fermi 
gas in a bag, by the same factor as the bag constant: 
Poi&c) = (1 — 2a c /7r) po(0). Since the stellar mass scales 

— 1/2 

as p , this implies that the least upper bound on the 
mass of the star as a function of the QCD coupling con- 
stant is given for non-rotating strange stars by 

-1/2 



M„ 



,(Oc) = 1 - 



2a r 



M max (0) . 



(3) 
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Fig. 3. The maximum mass of static strange stars in the 
MIT bag model, as a function of the strange quark mass, 
for various values of the QCD coupling constant. M max 
is always above the dotted line, which corresponds to the 
upper bound on B (Section 5); and below the continuous 
line, which corresponds to the lower bound on the bag 
constant, B min = 58.9(l-2a c /7r) MeVfm" 3 (eq. [f|). Also 
shown are lines corresponding to other fixed values of the 
bag constant: B = 80, 70, 60, 50 , and 45MeVfm -3 . 



through first order in a c (Prakash et al., 1990). For a c = 
0.4, eqs. (2), (3) give a maximum strange star mass of 
2.29A/q, higher by 16% than the maximum mass which is 
obtained for a = 0. 

Stellar rotation at a frequency up to 700 Hz would 
increase the maximum mass by only a few percent. In 
Fig. 2, we present the maximum mass of a strange star as 
a function of the rotational frequency, / = 1/ P. This is not 
the frequency dependence of the mass of a single star — 
what is plotted is the termination point, at each frequency, 
of a sequence of stellar models varying in mass. Thus, 
the baryon number of the maximum-mass model varies 
with frequency. In fact, for all / > 0, the maximum-mass 
stars are supramassive — if spun down at constant baryon 
number to / = 0, they would become unstable to collapse 
(compare the discussion of neutron-star models in Cook et 
al., 1994). We obtained the results presented in Fig. 2 with 
an accurate code based on spectral methods, developed by, 
and described in, Gourgoulhon et. al. (1999). 



4. The maximum mass of non-rotating strange 
stars as a function of m s and a c 

To determine the maximum mass of such stars for spe- 
cific values of model parameters (B,a c ,m s ), the e.o.s. of 
matter composed of interacting, massive quarks must be 
considered, and it can only be determined numerically 
The maximum stellar mass then following from the lower 
bound on B implicit in cq. (|^) is exhibited in Fig. 3 (con- 
tinuous lines) as a function of the strange-quark mass and 
of the QCD coupling constant. 

From Fig. 3 we can conclude for example, that if a 
strange star of two solar masses (2M Q ) were identified, 
and if a c were less than 0.4, then m s < 200 MeV (with 
both quantities rcnormalized at 313 MeV). 

If the density at zero pressure of quark matter is indeed 
close to its lowest bound given by the stability limit of 
eq. (0), then for a c > it can attain a value lower than 
4 x lCr 4 g cm~ 3 and the mass of even a non-rotating strange 
star may be sufficiently high to allow an orbital frequency 
of 1.06 kHz in the marginally stable orbit. It would appear 
that in the model of interacting quarks considered here, 
and in contrast to models in which po > 4.2 x 10 14 gcm~ 3 , 
the possible presence of moderately rotating strange stars 
in LMXBs could be com patible with the ke plerian model 
of kHz QPOs— compare |Bulik et al. (1999a| ) . 



5. The maximum value of B 

There is another stability constraint limiting the param- 
eters of quark matter, if such a substance is indeed the 
ground state of matter. The chemical potential at zero 
pressure of quark matter in three flavours and electrons 
in beta equilibrium, should be less than the rest energy per 
nucleon of the iron nucleus ^ d ' s < /i( 56 Fe) = 930.4 MeV 
(Farhi and Jaffe 1984). For massless quarks this corre- 
sponds to B/(l - 2a c /7r) < 91.49 MeV, but in general the 
constraint depends on the strange quark mass, and the 
maximum stellar mass corresponding to this upper bound 
on B is exhibit ed in each panel of Fi g. 3 as a dotted line 
[for details see Prakash et al. (199C ) and Zdunik et al. 
(2000)]. Note that the upper and lower constraints on B, 
when taken together, exclude high masses of the strange 
quark for low values of a c , if strange quark gas is to be 
the stable form of matter. 

6. Summary 

We have investigated the question whether moderately ro- 
tating strange stars with masses somewhat higher than 
2M© are allowed by relativistic equations of stellar equi- 
librium, and found that the answer could, in principle, be 
positive within the MIT bag model of beta-stable quark 
matter. An extension of this discussion to other models of 
strange matter will be given elsewhere. 

The physical constraints on the bag constant follow- 
ing from the basic hypothesis of stability of self-bound 



Zdunik, et. al.: On the mass of bag- model stars 



5 



quark matter (eq. M) allow B to be so small, that the 
corresponding mass of the star could be as high as 2.5M Q 
(eq. [3], Fig. 3). However, the strange star mass cannot be 
higher than 2.6M Q , even for stars of rotational period as 
short as 1.6 ms. 

The perturbative approach used here is sensible only 
if the value of the QCD coupling constant a c is small — we 
used < a c < 0.6. In this range, the direct dependence 
of M max on the coupling is practically negligible (Fig. 1), 
but the window of allowed values of B does depend on a c 
(Fig. 3) . The actual value of B is subject to a very large un- 
certainty. Fits to the hadronic mass spectrum (DeGrand 
et al, 1975) gave B = 59 MeV fm" 3 . 

The stellar mass decreases with increasing mass, m s , 
of the strange quark, and is lower by 10% to 20% than the 
one for massless-quark matter for the typical range consid- 
ered in the literature, 150 < m s c 2 /MeV < 300 (Madsen 
1999). Finally, the maximum stellar mass for strange stars 
rotating with a frequency up to 700 Hz, is larger than the 
one for non-rotating stars by less than 4% (Fig. 2). 

The physical reason for which large values of mass are 
obtained for a low value of the bag constant, is that the 
latter is proportional to the density of quark matter at 
zero pressure, po, while the maximum mass of the star is 

— 1/2 

proportional to p g . Thus, within the MIT bag model of 
quark matter, strange stars with mass M > 2.0M Q must 
have surface densities p < pi = 4.2 x 10 14 g/cm 3 (unless 
rotating with periods < 1.5 ms). 

This research was supported in part by KBN grants 
2P03D00418, 2P03D02117, 2P03D04013. The numerical 
calculations have been performed in part on DARC com- 
puters purchased thanks to a special grant from the SPM 
and SDU departments of CNRS. 
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